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We examine the importance of second order corrections to linearized cosmological perturbation 
theory in an inflationary background, taken to be a spatially flat FRW spacetime. The full second 
order problem is solved in the sense that we evaluate the effect of the superhorizon second order 
corrections on the inhomogeneous and homogeneous modes of the linearized flucuations. These 
' second order corrections enter in the form of a cumulative contribution from all of their Fourier 

, modes. In order to quantify their physical significance we study their effective equation of state 

■ by looking at the perturbed energy density and isotropic pressure to second order. We define the 

energy density (isotropic pressure) in terms of the (averaged) eigenvalues associated with timelike 
(spacelike) eigenvectors of a total stress energy for the metric and matter fluctuations. Our work 
suggests that that for many parameters of slow-roll inflation, the second order contributions to these 
energy density and pressures may dominate over the flrst order effects for the case of super-Hubble 
evolution. These results hold in our choice of first and second order coordinate conditions however 
we also argue that other 'reasonable' coordinate conditions do not alter the relative importance of 
the second order terms. We find that these second order contributions approximately take the form 
of a cosmological constant in this coordinate gauge, as found by others using effective methods. 
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INTRODUCTION 



|0 ' The major success of inflationary cosmology is in simultaneously offering an explanation for the homo- 

I geneity of the universe along with a mechanism that explains its inhomogeneity. The mechanism generating 

O . inhomogeneities involves quantum fluctuations in the fields that represent the dominant form of stress-energy 

during the inflationary epoch, which is a postulated 'potential' dominated era of the early universe (which 
$H . explains the homogeneity). The dynamics of the transition from the inflationary era to our current 'kinetic' 

or 'rest mass' dominated era are crucial in the formation of inhomogeneities, and in particular fix how much 
these quantum fluctuations will be amplified during the transition. In order to fully quantify the effect of 
these fluctuations on the large-scale geometry of the universe we can in principle use the field equations of 
. ^ general relativity to tell us exactly how the coupled matter and metric fluctuations clasically behave. Of 

. course, in practice, we cannot solve the full fleld equations for this scenario so we solve simplifled approxi- 

mations of the field equations within the framework of cosmological perturbation theory. At linear order the 
quantized cosmological perturbation theory has emerged as the primary tool to investigate the behaviour of 
fluctuations in the inflationary era, where e.g. it predicts an approximately scale-free spectrum of density 
fluctuations (see 01 for a comprehensive review). 

However, the beguiling mystery of the cosmological constant and dark matter/energy problems only deep- 
ens with increasingly accurate observations. There has been renewed interest in the past few years in the 
effect of higher order corrections to the linearized Einstein equations on both early and late-time physics 
in inflation. Suggestions have been made that superhorizon higher order corrections to linearized theory, 
on superhorizon scales (i.e. larger than the Hubble radius), can take the form of a negative cosmological 
constant. This could produce a dynamical relaxation mechanism for the bare cosmological constant (starting 
from 0, to 1^ most recently, and many references therein). The measurability or physical reality of 
these superhorizon backreaction effects has been a contentious issue (see e.g. jH, 0| and reference therein), 
and many questions remain regarding the link between local subhorizon physics and these superhorizon 
backreactions. 

In this paper we do not add to the discussion of measurability, but focus rather on explicitly evaluating 
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the second order corrections to the right hand side of the homogeneous Einstein equations in perturbation 
theory. At second order in perturbation theory one might even expect that, since the effect of second order 
contributions are cumulative over aU wavenumbers, their relative amphtude may become comparable to that 
of first order. Furthermore, if one thinks about solving the perturbed Einstein constraint equations for the 
matter fluctuations and putting these solutions back into the perturbed evolution equations, it is not hard 
to see that some of the second order corrections could in fact be divided by a so-called slow-roll parameter. 
This only adds to the worry that the second order terms could plausibly dominate for a 'slow enough' roll in 
the background. It is clear that in the limit as the slow-roll parameter goes to zero, so that the background 
universe tends to de Sitter, the first order corrections go to zero and the second order fluctuations dominate 
in their effect on the gravitational field. At what values of the slow-roll parameter do the second order 
perturbations dominate over the linear ones? The most radical possibility is that the slow-roll conditions are 
precisely the conditions that the second order perturbations dominate. Such questions appear to be behind 
some of the concerns raised by L. Grischuk in about the consistency of linearized perturbation theory in 
inflation. 

The technical complications of sorting out the second order gauge issues and other nonlinear effects such 
tensor perturbations seeding scalar perturbations are many, but rendered tractable with the aid of packages 
such as GRTensor for Maple llQj . In this paper we calculate the cumulative second order contributions to the 
homogeneous energy density and pressure. We do not address questions of the ultraviolet (short wavelength) 
regularization of the fluctuations but focus on the superhorizon fluctuations. Indeed, we pay special attention 
to the case where one considers the cumulative effect of Hubble sized to nearly homogeneous contributions 
on the homogeneous mode. 

In past work effective approximation methods have been used to evaluate such contributions. One popular 
method characterizes the backreactions in terms of an effective energy-momentum tensor Tat- In this method 
there are two contributions to Tat- the quadratic matter energy momentum tensor and the contribution of 
the first order gravity perturbations. Using early work by Brill, Hartle and Isaacson (0, 0) among others, 
the Einstein equations are expanded to second order assuming the linearized equations hold (so that they 
drop out). Then the remaining terms are spatially averaged with respect to the non-dynamical background 
metric and the resulting equations are interpreted as equations for a new homogeneous metric gat which 
include the effects of quadratic linear perturbations: 

Gabidab) = K{Tab+ < Tab >), (1) 

where < Tab > is the spatially averaged 'backreaction' stress-energy defined by 

Tab - Tab[iSgcd?,mf]~-Gab[{S9cd)% (2) 

K 

Here, k = SttG in units where c — 1, 5" indicates the n-th order perturbation of the object it acts on (as 
explained more precisely in the next section). 

It should be noted that, by construction, the zeroth order equations are not obeyed in this formulation. 
In this effective scenario one is solving for a new isotropic background metric which obeys equation (1), and 
this can in itself raise difficult questions of consistency if one is interested in backreactions on inhomogeneous 
modes (but is ok if one looks at just the homogeneous mode). This can be seen by considering the first 
variation of the Hamiltonian action for a gravitational system, namely 

SH = J {6N'Ha_ + SNiTC + N6n± + N^n')d^x, (3) 

and assuming that the background constraints do not hold, i.e. Ti^ 7^ 0, 7^ 0. If one has homogeneous 
variations then these linear terms will vanish anyway, whether or not the background equations of motion 
are satisfied. This is so because one can do a by-parts integration in the first and second terms above whose 
result will be boundary terms which can then be set to zero under reasonable assumptions 2] . 

Furthermore, it should also be noted that in their effective approach the full second order Einstein equations 
are not solved, nor are second order coordinate transformations considered. This latter fact can be of 
considerable concern when interpreting the significance of higher order effects . In this sense the effective 
approach does not appear, to us, to be able to convincingly evaluate the higher-order corrections to Einstein's 
equations, simply because it never actually considers them in the context of higher-order perturbation theory. 

Nevertheless, in and |3| this effective approach is used to evaluate the dominant long-wavelength 
contributions to Tab- Defining the energy density and pressure at second order by (using (+, — , — , — ) as the 
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signature) p =< t\ > and S^p = —(1/3) < t\ >, they find that these contributions have the effective 
equation of state S'^p « —S^p, with 6^p < 0. This corresponds to the equation of state of a negative 
cosmological constant. They also find that 6^p grows with time, partially because as inflation proceeds more 
and more length scales exceed the Hubble scale and contribute to S^p. These two results combined suggest 
their main claim, which is that the backreactions effectively create a negative, and growing, cosmological 
constant which can reduce the actual cosmological constant in the large. Locally one might expect the 
situation to be different. The first order modes should, locally, look like simple coordinate transformations 
of the homogeneous solutions, and their effect on higher order metric and matter fluctuations to be again 
that of higher order coordinate transformations. However, as we stated earlier, we will not focus on these 
difhcult issues of interpretation in this paper. 

Instead, we focus on evaluating the higher order contributions to the background equation of state, and in 
particular calculate the quantity S^p + S^p and its dipsersion. We follow a procedure of consistently (though 
probably not convergently) expanding the Einstein equations to second order and solving a subset of them 
assuming the zeroth and linear order equations hold. We do this about a flat FRW spacetime in which the 
dominant gravitating matter is a slowly-rolling, minimally coupled, scalar field 0, and we only study the 
effect of fluctuations on spatial scales exceeding the Hubble radius. In order to define the second order energy 
density and isotropic pressure, we give an invariant definition of 6^ p and S^p in terms of the eigenvalues of the 
stress energy tensor. Due to the mixing of tensor and scalar waves, these fluctuations will not only arise from 
second order scalar modes but also from quadratic combinations of scalar-scalar and tensor-tensor modes 

at second order. Wc find that in general 6^p + S^p ^ and 6^p < 0, but that ^ ^-t^ ^ does become small. 
Perhaps surprisingly, we also find that the relative amplitude of the second order dipsersion < (6^ pY > 
dominates over its linear counterpart < {SpY > for a wide range of slow-roll parameters in the background. 



In linear cosmological perturbation theory [3| the tensor, vector, and scalar metric modes all decouple, 
which justifies examining only one class of modes at a time to make the calculations simpler. Furthermore, 
depending on the initial conditions one takes for the linear metric modes at the beginning of inflation (a 
contentious issue), the scalar modes are typically more important than the gravitational modes at the end 
of inflation. Given this, the fact that scalar modes directly lead to energy density fluctuations required for 
structure formation, and the fact that for longwavelength perturbations gravity wave terms are typically 
suppressed by factors of fc^/a^, 'cosmological perturbations' have become synonymous with scalar pertur- 
bations. At second order, however, it is well known that pure second order modes have products of tensor, 
vector, and scalar modes as sources. For example, linear tensor fluctuations can induce second order scalar 
modes via the second order field equations. In order to describe this it becomes essential to include all three 
classes of modes (scalar, vector, and tensor) as sources for the higher order gravitational radiation. 

For fluctuations about a spatially flat FRW background in comoving coordinates {t, x) the perturbed 
metric may be written down as 

ds^ = -(1 + €A{t,x) + e'^A{t,x))dt'^ + 2{eBi{t,x) + e^Bi{t,x))dtdx^ + a'{t){5ij + ehij{t,x) + e qij{t,x))dxUx\'i) 

where e is the strength of the linear perturbation, i,j,k,... denote purely spatial indices, and a{t) is the 
usual scale factor. Together with the matter perturbations of the scalar matter (f), defined by 



the perturbations are {A, Bi,hij,^) at linear order and {A,Bi,qij,J-) at second order. The backreactions 
or higher corrections are suitably integrated quadratic combinations of terms from the former set, and they 
affect the longwavelength part of the latter set of variables. In principle, their evolution is determined by 
substituting the above metric and matter perturbations into the Einstein equations and solving these at 
second order subject to the linearized and background equations, i.e. we solve 



II. 



LONG- WAVELENGTH SECOND ORDER PERTURBATION THEORY 



0(i,x) = (t){t) + e<^>{t,x) + e'^T{t,x), 



(5) 




(6) 



and also demand that the linearized and zeroth-order (background) field equations 

6Gab = K (STab) 



(7) 
(8) 
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hold. Here, Gab is the usual Einstein tensor and Tab = {<l^;a<P;b — 9ab (5^''^^;c + is the stress-energy 

for a minimally coupled scalar field. 

In the following we assume that the spatial dependence of the fluctuations is of the form e='=''=*^', and we 
express the longwavelength approximation by 

Here, H is the Hubble parameter that corresponds to the scale factor expansion a{t) t", a ^ 1 and we 
take the potential of the slowly-rolling background scalar field to be 

= A + /3<^, (10) 

for which the only non- trivial slow-roll condition is ^(-^)^ = ^ 1- One should note that H = —^jj, 
and we write ^/a. ^ —^= = \/ ■^^^ so that the deSitter limit corresponds to a ^ oo. Conservation of 



energy for a slowly-rolling scalar field in this potential, with initial value (j)o, requires that (j) take the form 



(11) 



= <l>o-^, 

and we shall only consider comoving times t such that < t <C . 

Therefore we have three small parameters in this problem: e (the strength of the matter and metric 

fiuctuations), (slow-roll parameter associated with our choice of infiaton potential), and the long- 
wavelength parameter. However, only two of these small parameters arc independent as the order of the 
metric and matter fiuctuations is a direct product of the physics of the slow-roll parameter, so the scale of e 
is in some sense dynamically set. Unless otherwise specified, all future references to the order of a quantity 
will refer to its order in e. 

The solutions to equations (6) and (7) will be invariant under under a class of diffeomorphisms which 
are themselves functions of e. By expanding out the diffeomorphisms order by order in e, one can apply 
first order, second order, etc., coordinate transformations to the solutions of the first order, second order, 
etc., solutions of (6) and (7). In particular, one can choose these transformations to not only simplify the 
form of the solutions and the equations themselves, but to also deduce what part of the metric and matter 
fiuctuations is physical and what part is just a coordinate effect. 



Linear and second order coordinate transformations; gauge fixing 



As compared to the general covariance of the full (infinite order) theory, which allows for arbitrary coor- 
dinate transformations, within the framework of second-order perturbation theory we consider only linear 
and second order infinitesimal parts of these coordinate transformations (called gauge transformations in 
cosmological perturbation theory). The individual perturbations of the metric or stress-energy components 
will change in some well-defined way under such gauge transformations, and there will in general exist (an 
infinite number) of combinations of fiuctuations which are invariant, to second order, under this restricted 
class of coordinate transformations. 

Indeed, one can write a general e-dependent coordinate transformation by 



= i"(.T,e), 

and define an associated linearized coordinate transformation by 

dx''{x,e) dX''{x{x,e),e) 



C = lim 



(12) 



(13) 



where X"{x{x,e),e) 



de djf 

x", Ve, and similarly we can define the second-order coordinate transformation by 



e^o de 



dx''{x,e) dX''{S:{x,e),e) 



de 



dx" 



(14) 
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Since the metric gab is a tensor it will in general transform by 

9abix,e) = ——^gcd{i{x,e),e), (15) 

one can show that under a linear coordinate transformation the linear metric fluctuations Sgab suffer the 
change (using equation (13)) 

lim^^^^^ = Sgab ^ Sgab + £c{9ab{e = 0))^ gab + £cgab, (16) 

where is the usual Lie derivative along the vector (^°'. Similarly using equations (13) and (14) and the 
notation of equation (16), it is not hard to show that 

S^gab = S'^g'ab + i^l + £x)9ab + 2£fSgab, (17) 

so that in particular one can see that the change in (gauge transformation of) the second order fluctuation 
S^gab depends on C", i.e. it depends on the linearized coordinate transformation as well as the second 
order transformation x"'- This implies, for example, that the changes in the second order stress energy S^Tab 
of equation (6) caused by the gauge transformation will depend on both and x": the first and second 
order gauge transformations. 

In standard cosmological perturbation theory one usually makes a particular choice of (^"^ to simplify the 
interpretation of the fluctuations, e.g. the longitudinal gauge 

C° = B-aE, (18) 
C = ~d'E, (19) 

so that Sgoi — and Sgij cx gij . It should be noted that longitudinal gauge only fixes the scalar part of 
the metric into diagonal form, and its simple form relies crucially on the form of the anisotropics of the 
perturbed stress energy. However for our problem, since we obviously cannot diagonalize the entire linear 
order metric because of the presence of TT gravity waves, we find it convenient to emulate the harmonic 
gauge from the full theory, namely 

1 



da[V-\9\9'"']^0 (20) 



To linear order this gauge choice corresponds to setting 

b:. = 0, (21) 

a. (h,^ „ ^-^S,,) = 0, (22) 



which fixes, to within trivial residual gauge freedoms, the linearized metric perturbation regardless of the form 
of the perturbed stress energy and is sometimes known as the Poisson gauge Clearly the longitudinal 
gauge for the scalar sector is a special case of (21) and (22), since it restricts only the potentials B, E of the 
metric fluctuations B i^ The primary physical advantage in using this generalization of the longitudinal 
gauge is that one can unambiguously transform to any gauge while easily keeping track of the residual 
freedoms, while the primary mathematical advantage is that they lead to a compact form for the perturbed 
hamiltonian and momentum constraints. One can always transform from this gauge to any other gauge 
since the transformations are algebraic in nature (as opposed to, say, the nonlocal integrals that take one to 
synchronous gauge) (see !T5| for more details). 



III. TOTAL ENERGY DENSITY AND PRESSURE AT SECOND ORDER 



How do the classical metric and matter fluctuations at second order influence the background equation of 
state, and in particular how does one assess the influence of the gravitational backreactions? Fortunately, 
within flnite order perturbation theory we can avoid the conceptual and technical problems involved in 
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defining local or quasi-local definitions of gravitational energy-momentum because, by definition, we have 
a preferred decomposition of the spacetime metric. This allows us to exclusively attribute 'energy' to the 
'dynamical part' {5gab,5^gab) of the metric as opposed to the 'background part' {(jab) simply because there 
exists a 'background derivative' (Va) against which to measure any such 'dynamics'. Thus we can define a 
'relative gravitational stress energy' Tab of the fluctuations with respect to the curved background, and in 
particular we can take combinations of Tab and the stress energy Tab to study the fluctuations in the pressure 
and energy density at second order. 

A general formalism to define conserved quantities and conservation laws with respect to curved back- 
ground spacetimes has already been developed by Katz, Bicak, and Lyndon Bell in 13]. The basic idea is 
to start from the Lagrangian 



where 2Ca6c(e) = 2{Tabc{e) - f ahc) = bgca{() + ^ cgba{() - ^agbc{e)) and Lm is the 'matter' Lagrangian. 
Note that C^^ = 0, so that in the background L = 0. From L one can build vector densities which are 
conserved in the sense that V^/'^ = 0. We refer the reader to ,13] for the details of their construction and 
simply quote the result 



C (24) 



where is the arbitrary, smooth, vector field and is the analogue of the Einstein pseudotensor (defined 
in terms of C^^ instead of L^^). The first term on the right hand side of (24), in square brackets, can be 
interpreted as the relative stress energy of the fiuctuations with respect to a given background if we expand 
it to the desired order in e. The last group of terms can be interpreted in terms of the relative helicity of the 
perturbations with respect to the background and we shall not consider these in this work. We denote 
the gravitational parts of the relative stress energy by Tab, which is only a tensor to second order in e. 



A. Eigenvalues of the total stress energy 

One often defines (the rotationally invariant but not boost invariant) energy density and isotropic pressure 
of a perfect fluid by —p = g^^Too and 3p = g^Tu. While this is relatively straightforward to interpret and 
implement in linearized perturbation theory, at second order one gets complications such as having to subtract 
off shear (offdiagonal) stresses ~ Spi5pj,i ^ j, from the diagonal isotropic contributions ^ {Spi)^. From this 
and other points of view it turns out to be extremely useful to consider the eigenvalues of the mixed- valence 
total stress energy of the fluctuations, i.e. to consider the eigenvalues of the tensor 

T% + STl + 5^T% + S^T%, (25) 

where the last term <5^r°j is the second order part of the relative gravitational stress energy. Tab, described 
above. Since it is of mixed valence, it transforms from coordinates a; to x as {dx / dx){dx / dx) and therefore 
has gauge-covariant eigenvalues \i associated to timelike and spacelike eigenvectors, which are calculated by 
solving the equation 

det {f% + 5T1 + 5^T% + 5^T% - X,5\) - (26) 

At zeroth order, since we are not in deSitter but in a slow-roll spacetime, the one timelike eigenvalue and 
three spacelike eigenvalues (associated with their respective eigenvectors) are different and are in that sense 
'sufficiently separated'; this is sufficient to guarantee that the perturbations obey this property of being well- 
separated as well (the three spatial eigenvalues are not, but we only are interested in their average value). 
We define the energy density as minus the eigenvalue of the timelike eigenvector and the cumulative isotropic 
pressure as the average of the distinct eigenvalues associated with their respective spacelike eigenvectors. We 
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also emphasize that the eigenvectors at second order will in general point in different directions, but that 
all that matters in our calculations is the averaged contribution obtained after quantum averaging over k. 
After such averaging the terms like 5p{k)5p{k') collapse to the diagonal terms {6p(k))'^. 

We can express these eigenvalues in terms of scalars formed from the stress tensor and powers thereof. For 
example, to linear order one may find the averaged eigenvalues 6p and J2i ^Pi by perturbing the expressions 



5a ab 
a ~ 



-p- 



E 



p 



Pi 

2p 
3 



where S^b = T^b 



-f^gab- To linear order, these relations are equivalent to 
5{Tl) = -5p+^5pi 



(27) 
(28) 

(29) 

(30) 



where p + p = • Substituting the explicit expressions in terms of the metric and matter fluctuations 

for the left hand side, one finally obtains the desired expressions for the energy density and cumulative 
isotropic pressure. To higher order this procedure becomes more complicated since, as we mentioned above, 
pressure contributions like Cl2iPi)^ contain both diagonal contributions like J2iPi ^^'^ offdiagonal shear 

contributions like "^.j^j SpiSpj. However the basic strategy is similar and in any case equivalent to solving 
the above determinant using scalars such as the trace, 'double trace', and determinant, as we now sketch 
out. 

For the second order case, the starting point is 

S\T'^a+r"J = -S'p + Y,^'pi (31) 

1 

S\ls-bS\) = {5pf + ir£ 6p\ + E ^^'^ - ^ E ^Pi^Pi + + P)i^^P + ^ E '^'f (32) 



where Sab — {Tab + Tab) — "* 9ab- Thc appearance of terms like X^^^j SpiSpj and related cross-terms 

complicates the isolation of the desired eigenvalues '^^S^pi. In order to eliminate such terms we consider 
the second order perturbation of the cube of the trace-free part of the total stress energy. Combined with 
equations (31) and (32), this will give us another equation and with it the possibilty of cancelling these shear 
terms in terms of some function of metric and matter fluctuations. The general expression for the cube is 



-s%s',s\ = p^-J2pi-pl2pi + J2 p^p'm +p^T.Pi + ^pJ2p^pj - '^p^p^p^ 



(33) 



which to second order is 
= (P + 



and which in turn has the right form to solve for the shear terms we are desiring to eliminate from expression 
(30). Substituting in the expression of 5^(— |/S'°(,S"^„S''', ) in terms of metric and matter fluctuations and solving 
for the four averaged eigenvalues 5p, 5'^p, 5pi, Ylii ^^Pi in terms of the mertic fluctuations, we flnally get 



4{p + p){S^p + i S^pi) + i (Ei Spif + 3(V) + 2(5pEi ^Pi 



^S^P + E^S^P^ 

2{p + p){5p+\Y.,5pi) = 5i) 
-5p + J2i^Pi = ST 



(34) 
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The simultaneous solutions to these two sets of 



where 9 = ^SyS\, 6 = -lS\S',S\,T = g''\Tab + Tab). 
coupled equations are the (averaged) eigenvalues one would find directly from the matrix represented by the 
total stress-energy, expression (25). They are 



5p 

Sp + Sp 
S^p 
S'^p + S^p 



8(P+P) 



Pi 



1 

P+P 



SO 



_ 1 
P+P 

6(p+p) 



§11 4_ 
16 ^ 



24(p+j5) 



16(p+j5)2 



S-'T 
4 



(35) 



1. Linear contributions to the energy density and pressure 

At linear order the two values Sp, ^Pi of the total stress energy ST'\, comprised of only the stress energy 
of matter, can easily be found. Assuming the longitudinal gauge-fixing, we find 

Sp = -^(a,-3i7)a>-^Al .3g. 

Sp + Sp = 2/3$ J ' ^ ' 

To linear order, only scalar modes can induce energy density and pressure fluctuations. In the longitudinal 
gauge fixing specified by equations (18), (19), the (constrained) equation of motion for the spatial diagonal 
metric perturbation ^ in the long-wavelength limit, assuming slow-roll, is simply 

{d'i + Hdt)m = 0, (37) 

whose nondecaying solution can be taken to be a nonzero constant [ = ip ). This is to be distinguished from 
the pure deSitter case where this constant is precisely zero, since there are no physical linear scalar modes in 
pure dS . The corresponding matter perturbation $ is easily found via the constraint equations (namely 
$ = {3H'^//3)ip ). Using this and the constraint equations from (7) to express the result in terms of ip, we 
find the dominant contributions to Sp and Sp are 

<5p « ^ (38) 

^P + ^P - ^ (54eLM/ - {QHt - l)esR) , (39) 

where (^)^ = clw, -§t = esfl- Although the right hand of the latter equation is in some sense small, it is 
not zero. Therefore the linearized contribution to the equation of state, though highly suppressed, will still 
depend on the details of the small parameters. This point may perhaps be more obvious if one considers how 
the time evolution equation for the scalar field fluctuations ( during slow-roll ) is modified by the inclusion 
of the metric fluctuations, namely 

V 9(,^g^-9^); I 21 e,^) - (3 ' ' 

This shows how the gravitational fluctuations, at the linearized level, effectively induce a negative effective 
mass for the fluctuating scalar field as well as modify the effective Hubble parameter - all in a way which 
depends on what values we take for esB^^LW- Since the longitudinal gauge (or equivalently, a longitudinal 
choice of gauge invariant variables) admits no residual linearized coordinate tranformations, these fiuctuations 
cannot be associated with coordinate modes and are hence physical. 

It is also important to note that the right hand side of (36) contributes at a given, fixed, wavenumber 
k to the linear fluctuations Sp, Sp. At higher order we generically expect that contributions to the nearly 
homogeneous second order modes of S^p and S'^p will be cumulative over a broad range of k of the linear 
modes. It is this enhanced, cumulative, contribution to the second order modes that can make the nonlinear 
contributions to the equation of state nontrivial. 
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2. Second order energy density and pressure perturbations at fixed k 

Considering for a moment the situation at second order with fixed k, the partially gauge fixed eigenvalues 
of the total stress-energy are, for the generic case with no gauge-fixing, 



32k 



3iP 
"32^ 



+ 2^^ + 3F^^^^ + 18^ 



-6^p + 2I3T — + , 

4k k 



(41) 

(42) 



where is the squared amplitude of the linear tensor fluctautions. As in the linear case, the matter 
fluctuations are related to the metric fluctuations via the second order constraints from (6). The behaviour 

of the second order scalar perturbations will be influenced by not only the constant scalar modes at linear 
order, but also the linear tensor-tensor terms. One would thus expect that the scalar modes at second order 
will become time dependent. In fact, if we pick the 'harmonic-analogue' gauge of Section I, so that at linear 
and second order the scalar sectors of the metric are diagonal, we see 



t 



t + 
+ 

X X + a 



(43) 



B - a^'^'^E + 3E{3B - 10aE)k^ + A{a^E + HE){2aE - B) + 2aE{A - |oB) + B{aB - 2A) - 2aEi) 



-E + 



S^)E-—k'E'-2'4,E+-a'E'- 
2 ^2 



-B'^ - 4aEB 
2 



]} 



(44) 



then we obtain, using the spatial second order scalar field equations in this new coordinate system. 



A ■ 
B = 0, 

(2)^; = 0, 



-Q+-[{h+f + h+h- + ih. 



TT-TT sector 



^ V ' 

Scalar-Scalar sector 



(45) 

(46) 
(47) 



where qij = Qdij + {didi — Sij^Y^^'E and h+,h- € 3? denote the two TT independent degrees of freedom 
of hij. Note that the second-order lapse A contains contributions from the TT-TT gravitational wave 
contributions at linear order. In effect this is the longitudinal gauge at second order. It can easily be shown 
that it admits no residual scalar coordinate freedoms to second order. Within this gauge, under the slow-roll 
and longwavelength approximations and using the constraints to express everything in terms of the metric 
fluctuations, the equations of motion for Q are 



{df + Hdt)Qit) = 



2AHt - 162- 







(k/sO. 





70fc2 
"3^ 



9H 



(48) 



where = ((/i+)^ -I- h^h^ + (/i-)^) G SR. The growing solution for Q{t) is, for a fixed k mode of ip, C, and 
assuming that a ~ aot°' with a » 1, is 



Qit) 



24a2ln(t)- 



81fc^2 

K/32 



-2a 



(49) 



where D is a constant of integration deduced by the initial conditions of Q, implicitly set to zero by analogy 
with the linearized sector. This shows that the effect of a given linearized mode at some k on the nearly 
homgeneous mode of Q{t) is to make it time-dependent, as expected. 

When we compute S'^p and S^p in our exhaustive coordinate system we use the second and linear order 
constraints to eliminate the linearized matter fluctuation <I> and the second order matter fluctuation J^, so 
that the flnal expression is left only in terms of the metric perturbations. Following such a procedure, we 
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obtain the foUowing resuh for the dominant longwavelength contributions: 



'^(^^i/W + a.)Q(t)- ^^^^^ 



16k 



4k 



(50) 
(51) 



At first glance, it is immediately apparent that if we set C'^ — 0, Q{t) = (as Brandenberger et al do in 0| 
and references therein) then we obtain 5^p < Q but that S^p + is not in general zero because of terms 
that involve ratios of small parameters: 



6„/,2 



-54 



S^p + S^'p w +216 



fc2aV' 



a2t2«t4^2^2 ' 



(52) 
(53) 



However, we do find that, given = 0, Q{t) = 0, quadratic scalar-scalar backreaction contributions gen- 
erally do mimic that of a negative cosmological constant since {S'^p + 6'^p)/6^p ^ (k/aH)'^, which is set 
to zero in their analysis. Nevertheless, for sufficiently slow- roll it seems worrisome that for Hubble-scale 
size fluctuations the second order contributions S^p and 5^p can be comparable to 6p and Sp at a given k, 
since as they stand equations (52) and (53) show that the second order contributions carry 'extra factors' of 
slow- roll enhancement compared to equations (38) and (39). These extra factors come directly from solving 
the constraints for the matter variables [l^. 

However, it is clearly inconsistent to set Q{t) to zero in our approach since it is also of second order 
and in fact contains terms proportional to tp'^. Such a tactic also violates the second order field equations, 
which take into account that the metric fluctuation Q{t) at a given scale will receive contributions from all 
of the fourier linear modes under consideration (in our case, as we explain below, from the horizon to an 
approxiamtely homogeneous cutoff). Furthermore since the purely second order contributions to the energy 
density and pressure at these superhorizon scales are not constants but can evolve in comoving time, the 
contributions may change significantly during the slow-roll era to the end of inflation. 

In the following section we will examine the cumulative effect of the quadratic combinations of linear 
modes, averaged over superhubble scales, onto a given scale of the purely second order modes. For simplicity 
we will initially focus on the effect on the second order homogeneous mode. We strongly emphasize that 
the main gauge-fixing of the paper, as described by relations (21) and (22), is not well-defined in the 
strictly homogeneous limit fc = 0. An appropriate gauge transformation must be made to sensibly take the 
homogeneous limit and make these manipulations, and we discuss this in the next section. 



3. IR (super-Hubble) contributions from the backreactions 

We now consider the cumulative contributions to the energy density and pressure at second order due 
to the superhorizon modes. Considering only super Hubble fiuctuations we know that the dominant linear 
modes are independent of time since the background equation of state, during slow-roll, is approximately 
time-independent. In terms of the Fourier-decomposed ipk, the quantum fiuctuations (which we take to 
be Gaussian) during this era depend on k in such a way that the fluctuations per decade are a constant. 
Although strictly speaking cosmological fluctuations are quantized in terms of a reduced variable such as 
e.g. the Mukhanov-Sasaki (MS) variable v — a($ — -^jp-tp), one can always use the linearized constraints to 
simply relate (in the longwavelength limit) v and -0 up to time-dependent factors (see for more details). 
Indeed, the spatial two-point correlation function of ijj is (after an angular integration) 



< 0\4>{t,x)i>{t,x + r)\0 >=< Tp^ 



iH 



dk sin{kr) 
k kr 



r 1.3 



47r2 



(54) 



where ip(t, x) is the quantum operator associated with ?/), expanded in the classical basis of plane waves. The 
Fourier transform of the two-point function is the power spectrum, and completely characterizes Gaussian 
fluctuations in the sense that all higher correlation functions can be expressed in terms of it. Here, |0 > is 



11 



the vacuum chosen so that the modes of the reduced MS variable i^k obey Vkito) ~ /c"^/^, Vk ik}^"^ at some 
initial time Iq, which in turns implies a set of more complicated conditions on ^l^k which are not illuminating 
at this stage ( see Q for more details ). The metric fluctuations < -0^ > at the horizon scale are related 
to the density contrast fluctuations < {Sp/pY > by equations (38) and (39), and one can easily show that 
< {5pY >= 36(iJ'*/K^) < tjj'^ >. Using all of this, it is relatively straightforward to show that 



4(27r0)2 4esfl \2tt 

( see for more details ). We take this to hold to some almost homogeneous scale, say k = kmin « o,H, 
or in other words we cut off the infrared divergence of the linearized fluctuations at some scale k = kmin , so 
that for A; ~ 0, [i/'feP ~ 0. The factors of 4 on the right hand side come from using equations (38) and (39) to 
relate the fluctuations in the density contrast to the fluctuations in ijj. Again, this is equivalent to the usual 
statement made about the power spectrum in terms of the reduced Mukhanov-Sasaki variable v (namely 
that k^\i>k\^ ~ esRm'^ ' ^'^^ Pl^^*^'^ mass). It is worthwhile to notice that the corresponding result for the 
tensor amplitudes will not be enhanced by a slow-roll factor, so we will ignore them in what follows. 

Using equation (55) we can average over the quadratic combinations of linear fluctuations to solve the 
constrained equation of motion for the diagonal second order homogeneous metric fluctuation Qq (t) , equation 
(48). In other words, our strategy is to use the first and second order constraints to substitute for the matter 
fluctuations which appear in the Einstein equation for Qo(0; ^^'^ then solve for Qo{t) to calculate S^pand 

Indeed, keeping only the dominant terms and substituting for the matter terms, 

{d^ + Hdt)Qo{t) = 47r / ^^'>S{k)k^dk, (56) 

where ^'^'^ S{k) = 2AHt — 162 ^2*^2 (j^^ iV'fcP and kmin << aH. Carrying out the integral over k and 
solving equation (56), it is relatively straightforward to show that the dominant solution takes the form 

/272 \ 36iVKi72 aH ^ 

where N — J Hdt — Q:ln(t) is the number of e- foldings, and a — Ht >> 1 as described in Section II. Using 
this result and performing similar integrations for the remaining terms in equations (50) and (51), we finally 
obtain expressions for the contributions to the homogeneous mode of p and S^p 

S^PiB. {H + d„)Qo{t)~ / _-47r|^fe|2fc2dfc 



27^ 27, , , aH 



^3a(— - 6iV) + 36(7V - — ) ln( 



10 ' ' 2 ' ' fc„ 



(58) 



S'piR + S'piR w / {ATik^\'^k?)dk 

We can immediately compare the magnitude of the leading homogeneous backreaction term in, say, S^pm/p, 

to the root mean square of the density contrast < {5p/ p)^ > during inflation by using equation (55). One 
can demand a consistency condition for linearized theory, namely that the second order contributions be 

2 .[777:^ 



subdominant compared to that of the linearized sector, i.e. S pm < \/<{Sp) >. If we use the expressions 



(55), (57) and (58) above, this demand is crudely equivalent to the condition that, for N ^ 70, 



esR> {AnH^f'^N^'\ 



(60) 
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i.e. the usual slow-roll condition will in general be violated if [kH'^Y^^N^^^ > 1. Thus inequality (60) 
suggests, within our scheme, that the breakdown of the linearized approximation occurs when one assumes 
the slow-roll condition for the background spacetimc. Furthermore it is apparent that, although the right 
hand side of (57) is not zero, the form of these dominant contributions is approximately that of a cosmological 
constant since {S^p + 6'^p)/6'^p ~ -l/\n(^^fff-) ^ 0. 

As wc alluded to above, in the current gauge fixing we have chosen the homogeneous limit of the fluctu- 
ations is not well defined since, by equation (55) and the linear order equation ip = —A, the lapse A will 
diverge as A; — > 0. However, our results are valid even if one makes a gauge transformation which renders the 
superhorizon fluctuations well defined in the homogeneous limit. Indeed, after making such a gauge trans- 
formation, hen the homogeneous limit is taken and the total lapse goes to the value 1 while the offdiagonal 
terms go to 0, the central feature of 'slow-roll' enhancement remains and the above arguments still apply. 
Let us consider how this works explicitly by looking at the linear case. Choosing the infinitesimal gauge 
transformations to be = {T{t, k), kL{t, k),kL{t, k), kL{t, k)), we first note that Tk{t) must be, by equation 
(55) and the fact that Sq'qq = dgoo — 2T, so that 



ensures that \Ak\ for k ^ 0. Given this divergent gauge transformation of the comoving time, we 
can ask for what L]f{t) we can can ensure the rest of the perturbations that appear in the metric will be 
well-defined for fc — > 0. Since 

£c90i = -Ti + a''L,i, (62) 

in position space, then choosing Lk{t) = — / \/ Tsr'^^ ^^^^ render the shift B induced by equation (60) 

to zero. The off-diagonal terms E induced in turn by this choice of L will be well posed in the homogeneous 
limit since k'^E is what appears in the metric, and this will decay as Vk. However, the diagonal spatial metric 
terms will in general receive a contribution of 2Ha^T, which will get large as fc ^ 0. Since we implicitly use 
an IR cutoff of kmin, beyond which equation (54) will not hold and may in fact be replaced by a relation 
which does not diverge with k, these large contributions can be considered in some sense regulated. 

To second order the argument is very similar, only more tedious. Given the above choice for the linear 
gauge-fixing as fc ^ 0, wc must pick a second order x" = {^'^^T{t, k), 0, 0, 0) such that the second order shift 
and offdiagonal spatial terms go to zero in the limit fc ^ 0. Since 

S^g'oi = S^goi + SgoohT + Sguk'L + Sgf^'h^L + Sg,oif - 2k^L) + TdoSgoi - k^^^T + ^^^S 

where (2)5 = a?'kL{f - 12k'^L) - kT{At - AaaL - Qk'^L - a^'t), one can show that the choice 



(2)T = - 

k 



^^"i S + AkT + SipkL 



(63) 



will yield S'^c/'q^ — 0. The offdiagonal, second order, spatial terms induced by the above transformation are 

2k^T^ ^4ak^{2aL + aL)T+18a'^k^L^, (64) 

which at worst diverge as k~^. Since only the expression k'^{^'^^E) appears in the metric, these offdiago- 
nal terms in the metric smoothly go to zero in the homogeneous limit. Once again the spatial diagonal 
contributions to Q will grow large ask ^ Q. 

Now we can finally address the total effect of all these transformations on the quantities of interest, 5'^p, 5"^ p. 
To second order, for example, 5'^p will be Lie-dragged along C,"' and according to the tranformation 

= 5'^p+{£'l + £^)p + 2£(^6p (65) 



= 5y + 2{TdoSp - Sk^LSp) + ^^>Tp + [T^p + TTp - pk^LT 

and similarly for the averaged pressures 5'^p. Choosing the above expressions for T, L, (^^T, we can see that 
the total effect on S'^p and S'^p will be to introduce terms that diverge like k but are suppressed by factors of 
/?. Such terms will be subdominant compared to the terms already present in the original gauge, or in other 
words, the gauge transformation which renders the metric diagonal in the homogeneous limit does not undo 
the dominant contributions to the energy density and pressure at second order. 
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FIG. 1: All linearized modes in the shaded region, spanning from k = aH to k — kmin « clH in spatial scale, are 
taken to seed second order modes at a particular value of k. The thick black boundary indicates the Hubble scale, 
where the amplitudes of the linearized fluctuations freeze out during slow-roll. We compare the amplitudes of the 
linearized fluctuations at the thick black boundary to that of the second order fluctuations at = fc and k = (and 
ignore the influence of suitably renormalized subhorizon (k > aH) modes at second order) during slow-roll. 



B. Comment on backreaction on inhomogeneous second order modes 

An important loophole in the above analysis resides in the fact that we compare the horizon scale (frozen) 
amplitudes of the linearized fluctuations to that of the homogeneous sector of the second order fluctuations. 
It is far from clear that this is an acceptable comparison to make, not least because we are directly cutting off 
the divergence of the linearized fluctuations by imposing an IR cutoff at fc = kmin << o,H and then comparing 
the amplitude of the second order fluctuations with that of the linearized fluctuations well beyond the cutolf. 
A valid criticism of this result is thus that it would be natural for the second order perturbations to dominate 
at the homogeneous scale simply because we have cut off the linearized fluctuations long before comparing 
their amplitude to that at second order, or in other words the spatial dependence induced by evaluating the 
(quantum averaged) second order amplitude at some kmin < k = k << aH may alter the conclusions of 
inequality (60). See Figure 1. 

However, we find that when we indeed compare the quantity J < (^^~|^) > some scale k such that 



k << aH to the horizon-scale amplitude of the linearized fluctuations J < (^^^ 



>, we once 

k=aH 

again find that the amplitude of the second order modes may still dominate over that of the linear modes 
assuming slow-roll. As show in detail in the Appendices, we find that we may write 

< '5Vfl.(fc)'5Vi?(fc) > « < / / ( ) ^ik'-k)^k'^(k"-k)^k"d^k'd^k'' > 

+ <J (^^^ CQk'-kCQk'd^'k' > 

+ < [ [ ( — ^^]cQ{t',k';k)ij^k"-k)^k"d^k'd^k^'>, (66) 
where H = (9o + ) 7 find that this in turn has the form 

< 8^vinm''v\nk^) > ~ (—) 4^ (^1"' + " (^1 1^(7) + Ci ln(a))) , (67) 



14 



where 



2657205 „„„ 

-1594323 , , 

2657205 , , 

^ -68^ " 



and where 



7-1 + ^ (72) 

/I'm?,??, 




Using all the above we can now finally see that the demand that 
as described in Figure 1, is equivalent to 

esR>l{KH^)HA,N)-- (73) 

Given that ^ 70, Ai ^ 200, (73) implies e.g. that for kH"^ ~ 1 ^ ^ mpianck^ the slow-roll condition 
must be violated. We see that the amplitude of the second order flucuations in our quantity dominates 
the corresponding linearized amplitude if the background spacetime is rolling slowly enough and kH'^ is large 
enough, as it is in many models. Note also the appearance of the number of e-foldings iV, which indicates, 
as previously shown, that this effect is a cumulative effect which depends on the growth of the phase space 
of superhorizon modes. All in all, the violation of this inequality in such models casts doubt on the viability 
of the linearized approximation to those slowly rolling spacctimes. 

As we imply in the caption to Figure 1 above, this calculation would ultimately contain formally divergent 
subhorizon contributions to i5^p, 6^p and it is this real phsyical effect of the coupling between subhorizon 
(UV) modes with superHubble (IR) modes that will reveal the observable importance of backreaction effects 
for local observers. The points we are making here are 1) the superHubble contribution, in this reasonable 
gauge, docs take the approximate form of a cosmological constant contribution; and 2) assuming the slow- 
roll condition in this gauge, within a wide class of initial conditions, is equivalent to assuming that IR 
backreactions dominate the linear terms. Even if one uses the residual homogeneous gauge freedom, one 
cannot simulataneously gauge away 5"^ p and 5^ p. These results hold even if one considers backreaction onto 
some spatial scale k « kmin of the second order modes. 



IV. SUMMARY AND CONCLUSION 



To summarize, we study the perturbations of the inflationary slow-roll spacetime which are at second 
order in the metric and matter fluctuations. We follow a procedure of consistently (though probably not 
convergently) expanding the Einstein equations to second order and solving a subset of them assuming the 
zeroth and linear order equations hold. Specifically, we solve the first and second order Einstein constraint 
equations for the matter fluctuations, and in that way express our final expressions of and 5'^p in terms 
of the metric fluctuations. We then solve a second order Einstein to express all of our expressions in terms 
of quadratic combinations of linearized metric fluctuations. 

In order to isolate the physical degrees of freedom in the second order fluctuations we use a longitudinal 
gauge-fixing procedure at second order. Namely, we specify two independent infinitesimal, inhomogencous, 
coordinate transformations (gauges), one at linear order and one at second order, which admit no residual 
coordinate freedoms. Within this coordinate system we evaluate the fluctuations of two independent back- 
ground scalars formed from the stress energy and from them define the fluctuations, to second order, of 
the isotropic energy density 6^p and pressure 5^p. These fluctuations will not only arise from second order 
scalar modes but also from quadratic combinations of scalar-scalar and tensor-tensor modes at second order. 
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Whereas the nondecaying Hnear scalar and tensor modes are constants at these scales (the vector modes die 
away), the nondecaying second order modes are time dependent and this leads to time-dependent p and 
We calculate p and b^'p using the linearized and second order constraints and one second order scalar 
Einstein equation. 

Futhermore, given that we have three effective small parameters in this problem (a slow-roll parameter, the 
strength of the metric and matter fluctuations, and the longwavelength approximation (iJ^/A^ << 1), we find 
that our b"^ p and (5^p depend sensitively on the hierarchy of small parameters one assumes in the sense that 
ambiguous terms like {}? j H"^ c?){H^ j k0^) appear in the mode expansion of J^p, J^p. For the (incomplete) 
case of just scalar-scalar backreactions and no genuinely second order metric or matter fluctuations, we 
find that p < but that S^p + S^p 7^ in general. We find that the second order contributions to the 
energy density and pressure can, with the assumption of slow-roll, dominate over the second order linear 
contributions to the energy density and isotropic pressure given a broad range of slow-roll parameters. It 
is worth emphasizing that the divergence of the enegy density or pressure perturbations as the slow-roll 
parameter goes to zero is certainly not physical, as our expressions (say, equation (67)) would seem to 
suggest. Rather, if one thinks of a single set of Einstein constraints expanded order by order (as opposed 
to a tower of constraints) then it may well be that second order terms in these constraints simply start to 
dominate the linear ones. 

We conclude that when one truly goes to second order and solves some of the Einstein equations for the 
higher order classical fluctuations, they do approximately lead to a cosmological constant type of contribution 
in this gauge. Furthermore, it seems that these higher order corrections dominate the linear terms if slow- 
roll holds in the background, suggesting the breakdown of the linearized approximation to within a small 
window of slow- roll parameters. Some previous calculations of higer order super horizon effects (of which 
we are aware) have used a procedure which effectively takes the expectation value of the gauge-fixed metric 
before forming some sort of 'invariant' measure of the expansion with which to probe 'local' backreaction, 
i.e. gauge fixing before taking the expectation value. Such a procedure suffers from higher order gauge 
ambiguities, and at least for a model with massless, minimally coupled scalar with quartic self-interaction 
(no gravity) one can first form a desired 'invariant' and then take expectation values and gauge fix this result 
(see for example 0|, and also |l3| )■ Further investigations along this line will almost certainly prove useful, 
and one pay-off seems to be new ideas for observables in backreactions, such as recently described in [T6|. 

Finally, as we indicated in the Introduction, the procedure we use here does not say anything about what 
an observer would measure as the averaged cosmological constant in his own obserable patch of the universe. 
In other words, the suitably renormalized subhorizon contributions to S^p and 6^p will allow a probe of the 
possibly crucial physics of the coupling of subhorizon and superhorizon modes in inflation and possibly shed 
some light on what we even mean by local modifications to a cosmological constant. 
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APPENDIX A: FOUR-POINT CORRELATION FUNCTIONS 



The central result of this paper lies in comparing the second order quantity J < (^^-|^) > a.t some scale k 



such that kmin < k << aH to the horizon-scale amplitude of the linearized fiuctuations . /< 



>. 

k=aH 

The angled brackets indicates averaging over a Hadamard vacuum state |0 >. In this Appendix we supply 
details of the calculation of the four-point function which are inherent in the expressions for the dominant 

terms of < j >■ We first derive the expressions for the inhomogeneous case (k ^ 0) and then we take 

the homogeneous limit (fc = 0). The Appendix ends with a section which proves that second order scalar 
metric perturbations are non-trivial in pure de-Sitter (i.e., in a no- roll background), contrary to some claims 
made in the literature. 
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The case of the 2-point function is relatively simple to consider. Expanding ipk in terms of the creation 
and annihilation operators a, we write 

V'fc = t^fcOfe + 4,, (Al) 
where cok is equal to ( by analogy to the familiar Minkowski result u)k ~ (^TJi) ) 



k 



The first term is a normalization term that goes as 1/ Vv since corresponds to the volume measure of a 
comoving observer. This solution to equation (37) (with the spatial gradient term restored) is valid to within 
a couple of Hubble times of the horizon exit, and during this time the variation of H is negligible. Therefore 
it makes sense that, up to some phase factor that varies slowly on the Hubble timescale, the expression of cok 
has this simple form compared to the flat space result. 
Defining the 2-point function as 



/ il}(^k'-k)i'k'd^k' 
Jn,,, 



<il)^ >=< I ^^k'-k)i^k'd''k' > 
we find that 

<^2^ = </ (afe/_fcW(^,_g) +4,_fca>*g,_gP(afc'Wg, +wi4,)f^^fc"*' > (A2) 

= < ak'-kal,uj^p_^^u^l,(fk'> 

= S{-k) < j u^^,_^^u;tSk' >, 

where flk' indicates the superHubblc range of integration for k', namely k' E [kmin + k, aH], kmin << iH. In 
other words, the quantum average of the 2-point function is only non-trivial for homogeneous contributions, 
as suggested by the Poincare invariance of the linearized fluctuations. In the following we will assume that 
the spectrum of modes is discrete by imposing periodic boundary conditions {x x + L) on the spatially 
flat slicing of the background. In this manner the delta functions that appear in the commutation relations 
of ladder operators simply become Kronecker deltas and the volume normalizations are implicitly carried by 
the above definition of cuk- 

For the four-point function the situation is different. More specifically, we wish to compute the 'square 
of the two-point function', i.e. we do not want technically want the four-point function but rather the 
fluctuations in the operator ip"^ since we want to measure and compare the dispersions of p,p at first and 
second order. Indeed, denoting these fluctuations by < > we define (in the continuum limit) 



/ / ilJ(k'-k)ipk''^{k"-k)'^k"(fik'd^k" >, 



one may verify that 

<^^> = ^^k'-k)A'd^k' >^ (A3) 

+ < / / iJ(k'-k)'^{k"-k)d^k'd^k^' >< / ilJk'ipk"d^k'd^k^' > 

+ < [ [ i^(k'-k)^k"d''k'd'k^' >< j j xP(k"-k)^k-d:'k'd''k" >, 
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where similarly flh" indicates the supcrHubble range of integration for k", namely k" € [kmin + k,aH]. 
Assuming that A; ^ we find that, using equation (A2), 



< V > = < 
+ < 

which straightforwardly simplifies to 



ak'al„uj^,ujl„d?k'(fk/' >< / 



ak'-kalu_k^^j^,_k)^lp,_k)d''k'd?k!' > 



*^„(fk'(fk'' >< / ak''-kal,Lo^j;„_^^uj*^,(fk'(fk'' >, 

J 0,1./ JCli.// 



1 %n f H 



(A4) 



using equation (55) in the main text for k > kmin- It turns out that we shall also require the expressions 

< '4'(k' -k)''Pk''>P{k" -k)i'k" k'^ > and < '4'(k' -k)'4'k''4'(k" -k)'4'k" k"'^ k''^ > in what follows, so, using the above, we 
find (for A; ^ 0) 



< i^{k' -k)i'k'ip(k" -k)'4'k"k''^ > = 



+ k''\u.„\^d^k] ( [ 

+ ([ k" 



hk"-k)\'d'k'' 



k"\k»-k/d'k-^' 



\^Ck"-k/d'k^' 



■'(k"-k)\ 



(A5) 



and 



< i)(k'-k)'4'k'i^ik"-k)^k"k"'^k''^ > 



k"^\utSSk" 



k'\^k'-k)\'d'k' 



\u;j:„\'d'k^' 



-{k'-k)\ 



k'\k'-k/d'k' 



^^d^k' 



7 fc"Vfc"i'rf'fc''l 



k''\^,k,.k,?d'k' 



k"\i,k"-k)?d^k^' 



Hk"-k)\"d'k^' 



k" \^i^k"-k)?d'k" 



(A6) 



Once again taking equation (55) to give the amplitude of the linearized quantum fluctuations at horizon 
crossing and assuming this frozen amplitude for k' , k" « aH (up until k' , k" = kmin, where we cut it off) 
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we complete the calculation, using 



7r /'27r 



1 



1 9k / H 

lesR \2tt 



^0 (fc2 + fc' -2fcA;'cos(6l)) 

2 



3/2 



(i0 iim{6)d9 



1 



47r ( -(ln(fc'-fc)+ln(/c' + fc)) 



off 



and similarly 



/ 



i^f^l 4.(ln 



1 9/t /i? 



ln(l + 2^)), 



( — j 27r (fc'' + fc2(ln(fc' -k)+ ln(fc' + A;))) 



aH 

kmin + k 



(A7) 



(AS) 



1 9k / H 



AesR \2tt 

along with (remembering again that fc << aH) 

1 9k / H 



— 2tt[ {aHy + 2k' In 



aH 1 



ln(l + 2^) 



j k'^\uo^^,_j^^\''d^k' = \^ ^ {k'^ + '^k'^k'' + 2fc4(ln(fc' - k) + ln(fc' + fc))) 



aH 

kmin+k 



TT (aH)*+4F In 



ai/ 1 



1 9k f H\ 

to obtain finally (using the fact that the fc', k" ranges of integration are identical) 



fc . 2^^(' + '^^ 



< i'(k'-k)i'k''ip(k"-k 



)^k" > ^ T] (2(ln((T))^ - ^ ln(or)ln(7) + 1 (ln(7))= 



< ik{k'-k)'4>k'i>(k"-k)i'k"k''^ > » '7(o-H') 



1 + 2 



in 



2^7/ 2 vv^y 



(A9) 

(AlO) 
(All) 



< iJ(k' -k)tpk'ip(k" -k)'<Pk" k'^ k"^ > « '?(a-H") 



2 4 



a// 

In the above, we define the dimensionless factors (J,j,r] as 

aH 



l + l(W,) + 31n^)U('4Vfiln^ 



V7 



2 V7 



+ -?7 3 In— +lno-ln — 



V7 



V7 



1 + 



2fc 

2 



(A12) 

(A13) 
(A14) 

(A15) 



and once again we assume fcmm << aH. We retain powers of k/aH for now simply for generality. The 
long-wavelength approximation will kill off these terms later on in the calculation. 

Finally, one can take the homogeneous limit (fc — > 0) of expressions (A9)-(A11), bearing in mind that 
the squares of the two-point functions now contribute as shown by equation (A2). When we take the 
homogeneous limit of the above equations we find 



< V''' >fc=o ~ 3?7(lnCT) 
lim < ■!/'(fc'-/c)'0fe'V'(/c"-fe)V'fc"fc'^ > ~ 'dr]{aHfhia 

lim < V(fe'-fe)V'fe'^(fe"-fe)V'fe"'^'^fc"^ > ~ ^riiaHf 
which provides a coarse but useful check on the algebra to this stage. 



-l-lniT 



(A16) 
(A17) 

(A18) 
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APPENDIX B: SECOND ORDER EQUATIONS OF MOTION AND CORRELATION 

FUNCTIONS AT SECOND ORDER 

In order to compute the full quantity S'^ p{k)6'^ (k) we shall formally encounter two-pomt functions not 
only involving ipk, but also those involving the second order fluctuations Qk, e.g. 

< / CQk'-kCQk'd^'k' >, (Bl) 

where C = (da + H). Of course, by solving the second order equations of motion these sorts of expressions 
can be reduced to four-point functions involving only -0^. We now show how this reduction is accomplished. 
Using equation (56) and the fact that H — —H^^-f^, one can easily show via by-parts integration JJj that 

CQk ^ J '^^^S{t',k';k)d^k'dt' + esBH{t'fQ(t',k';k)dt'd^9, (B2) 

i.e. we compute the first integral of the reduced second order equation of motion. A tedious integration 
reveals that the leading terms of the latter integral over Q are suppressed by a factor of esn/oi compared to 
those of the first term. Therefore we ignore the latter terms and write 

CQk ~ f f ^'^^S{t',k';k)d^k'dt' (B3) 

which immediately leads us to the expression, again using equation (56), 

<CQ,> = < f j^^ Hit'f (^2AH{t')t' 162 ^^^,^,^2(t0^^iy ) ^i^'-'^)'^'^"^''^''^'' > (^4) 

Equation (B4) allows us to evaluate the relevant quantities which will appear in the expression for 
5^p{k)S'^p{k), such as (Bl) and < J^^^ J^^^^ CQ{t' ,k']k)'il)(k"-k)'^k"d^k'd^k'' >. Notice that CQ has units 
1/s, as makes sense since the metric fiuctuations are defined as dimensionless. Using the early results of 
Appendix A, we find that 

< CQu'-kCQk'd^'k' > = < j H^(^Am' -lQ2^j^^^ ,P(^k'-k)^k'd^k'dt' (B5) 

X 1* (^2AHt" - 162-^^^ ^P^k"-k)^k"d'k^'dt" > 

<[ [ CQit',k';k)^P^k''-k)A''d^k'd^J^' > = <[ H^2'iHt' -162^^^] [ ^(^k'-k)^k'd^k' dt' 



X / ^(k"-k)^k"d'k" >, (B6) 

where the appropriate time dependences of a{t) ^ aot°',H{t) ^ a/t are assumed above (as properly shown 
in equation (B4)). Expanding the above expressions (assuming once again that the ranges of integration for 
fc', k" are the same, as in Appendix A), we obtain 

< j CQk'-kCQk'd^k' > = J J H{t'fH{t"f(2Aft't" <ij^> dt'dt" (B7) 

' H{trH{t" f{2At') (^(^^) < V^fc'^ > dt'dt" 



ffHitrHm2u'')[-^)<,^k'->dt'^ 



t ft / 2 
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and similarly 



< 



/ / LQ{t\k']k)i}(k''-k)i^k''d^k'(fk'' > = 2A ( H'^im') <'4)'^ > dt' 
Jny Jay, J 



(B8) 



/ 



-162 / H 



, 1 H'^ 



< V^fc'^ > dt', 



where we define < ^p"^ >=< V'(fe'-fe)V'fe'V'(fe"-fe)V'fe" >,< tp'^k''^ >=< ^{k'-k)i^k'i'{k"-k)'4'k"k'^ > and < 
V'^fc" k" >=< ip(k' -k)''Pk''^l^{k" -k)4'k"k' k" > all integrated over k J k ^ as in the above sections. 

Inserting expressions (AlO) through (A12) (and Tising (A5) and (A6) to perform the temporal product 
integrations in (B7) [which amount to symmetrization in t' , t"] ) into (B7) and (B8) we obtain (for k 0) 



u/2 



< ^ CQy.kCQu'd^k' >= y* H{t'fH{t"f{2Aft't"f) (^2(ln(a))' - ^ ln(a) 111(7) + \ (ln( 

-/'/' 



(7))^ dt'dt" 



H{t'fH{t"f{24i') 



162 



-fj*H{t"fH{t')\2At")[^ 



fiiaHY 



1 + 2 



aH 

k 
ah 



V2V7y 2 vv7 



In 



V7 



In 



^/7 



i + i(ln(*) + 31n^) 



3fln4] +lnaln4 



V7 



(B9) 



and 



< 



^ CQ{t', k'; k)iP(k"-k)il'k"d'k'(fk'' >= 24 j' H^{Ht')ri (2{\n{a)f - | ln(<T) ln(7) + 1 (ln(7))^^ dt' 



1 + 2 1 



\aH) VVt/ 



In 



^).l.(^)).., ,B10, 



where we remind the reader that a, ^,7] are defined in definitions (A13)-(A15) and the "symbol denotes 
symmetrization in t' , t" . One can once again verify that the dimensions of all of the terms in equations (B9) 
and (BIO) are inverse seconds squared and inverse seconds respectively, as required. 

Once again we can immediately take the homogeneous limit directly from the above equations (again 

remembering to add in the contributions from the squares of the averages of the two point functions, as 
in Appendix A), or by using the relations (A16)-(A18) directly in (B7) and (B8). For example, the result 
for Imik^Q < Jq ^ Jq ^Q{t', k'; k)'tp(^k"-k)'4'k"d^k'd^k" > is, putting in the explicit time dependence of the 
background, 



lim < 

fe->0 



/ / 



CQ{t', k'; k)^^k„_k)A"d''k'd''k" > 



3(24) r^L^ 

"^^^^^ (27r)4 /^2^4 J ^,14 [^^ ^, 



,/a \ 2 



(Bll) 



-3(162)- 



(27r)4 K/32k2^4J ^ 



dt' , aoi'"a 
in 



and similarly for limfe_>o < J^^, CQk'-ki^Qk'd^k' >. 

The key issue is now one of simplifying (B9)-(B10) to extract their dominant parts (those 'most enhanced' 
by factors of the slow-roll parameter). It should be noted that this is complicated by the fact that the 
dependence on the slow-roll parameter is not just carried through factors of V' ^ f3, but also through a, 
which in the Introduction we defined to be proportional to 1/y/esR- Therefore we must, to be safe, carry 
out the time integrals before we make the approximation that a >> 1, esn « l,a/t ^ H G 3? epsecially 
since, jiist as differentiation in comoving time of the scale factor generally introduces factors of the slow-roll 
parameter in the numerator of a given expression, integration introduces factors in the denominator. This is 
the reason it is not immediately obvious the second term in equation (B2) is subdominant to the first term. 
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We begin in the simpler case of the homogeneous limit. In that case, one can for example show that the 
dominant terms for the following integrals are (for a » 1): 

* dt' / aot'^aV 2a^ 2a 



r dt!_ 
J t'^' 



and similarly for the other integrals required. Using these sorts of integrals in equations (Bll) for example, 
one can show that 

lim < / / £Q(t',fc';fe)V(fc"-fc)Vfc"d'feVfc"> = -3(162) ^^!^\^' f 3 [-a^ + a(-171n(a))] , (B14) 

and similarly for limfe^o < Jq CQk>^k^Qk'(Pk' >. 

For the full inhomogeneous problem we can make excellent use of the long-wavelength approximation to 
rewrite equations (B9) and (BIO) as 



< ^ rQfe,_fe CQyd'k' >= y* y" H{t'y'H{t"y\24ft't"fi (^2{ln{a)f - ^ ln{a) 111(7) + \ {Hl)f^ dt'dt" 

- f f HitrHit"n2M') (^) .ioHf (^) + im [j=)) dt'dt" (B15) 



t t-t 



i + l(ln(^) + 31nA) 



dt'dt" 



and 

< 



CQ{t', k'; k)'4>(u"-k)^k"d'k'd'k/' > K 24J' H\Ht')n (^2(ln(a))" - | ln(a) 111(7) + \ (ln(7))') dt' 

-16./'„»|;,(l.(j^)+ll„(^))..' ,B16, 
It turns out the dominant terms in these integrals are, for a >> 1 and esR « 1, 

< / CQu^-kCQu'd'k' > « -^i!^^^ [192a + 7040 In(cT)- 2112 ln(7)] (B17) 

</ / CQ{t',k';k)^(y,_k)'4>k"d''k'd'kj' > ^ .f,jfT 3 [32805q' - 6561q(-85 lii(a) + 51 111(7))] (B18) 
Ju., Jn.i, 4bJ47r esR 



where we note that k dependence comes in solely from the ln(7) terms and that the dimensions of equations 
(B17) and (B18) are respectively inverse seconds squared and inverse seconds, as required. 
We are now finally in a position to collect all of these results, namely, equations (B17)-(B18) and (AlO) 



(and (B15)-(B16), (A16)), to fully evaluate the quantity y < ^ ^^^" ^ > at some scale k such that kmin < 
k « aH (and also at fc = 0). The goal is to compare the magnitude of this term to the horizon-scale 



amplitude of the linearized fluctuations w < 



>, as shown in Appendix C. 

k=aH 



APPENDIX C: DERIVATION OF INEQUALITIES (60) AND (73) 

In order to derive inequalities (60) and (73) we must calculate the VEV amplitude of the second order 
fluctuations at the homogeneous scale. We shall derive the result for some scale kmin k « aH and then 



22 



take the homogeneous hmit to recover inequahty (60). Using the results and notation of Appendices A and 
B on the dominant (scalar-scalar) terms of equation (50), one may write the averaged square of the dominant 
contributions to the second order IR pressure contribution (or energy density, etc.) as 



< S^piR{k)S^j,{k) > « < / / 

/ (- 



KeSR 



V'(fe' -fe) V'fc' '>P{k"-k) V'fe" d k'drk" > 



+ < 

+ < 



flu, Jn 



CQk'-kCQk'Sk' > 
6H 54iJ2 



K KesR 



CQ{t', k'; k)i,(k"-k)i'k"(Pk'd^k" > (CI) 



Putting in the results from all of the above appendices we find that the dominant terms are of the form 



< S-'piR{k)S^j,{k) > « ( — 



where 



Ci 



^SR' 



{Aia^ + a {Bi \n{j) + Ci ln(a))) , 



2657205 
11567r2 
-1594323 

687r2 
2657205 
687r2 



(C2) 

(C3) 
(C4) 
(C5) 



It is important to note that the details of this calculation confirm that the naive guess ventured in the 

Introduction holds: there arc solutions of the diffcomorphism constraints of general relativity, at second 
order in perturbation theory about a slowly rolling background, which introduce a factor of 1/ csr into any 
expression of S'^pm and S'^pm, and these factors survive through quantum averaging. A second order gauge 
transformation cannot eliminate the presence of such slow-roll enhanced terms. 

The central result of this paper is in comparing the amplitudes of the second order fluctuations at kmin < 
k = k « (iH to the amplitudes linearized fluctuations at the horizon scale, k = aH. Using equation (C2) 
and equation (55) for this purpose, we see that 



< 1^ 
P 



> > 



k=aH 



k=k 



is equivalent to demanding that 



A2y/K.H 



y/eSRT^ 



> 



_KlP_ 



N 
esR 



N 

CSR 



{B, ln(7)+Ciln(a)), 



(C6) 



where a = N/esR, A2 = We can see that the spatial dependence is subdominant in the sense that it 
is multiplied by a lower power of a: in other words, the homogeneous fluctuations (mean) are not, in the 
slow-roll limit, altered by the finite corrections which occur when one evaluated at the finite spatial scale of 
k = k << aH. We also note that (C6) is positive definite, as it should be, since k << aH. Furthermore 
we note that the first and second order dispersions above have the same weighting in volume normalization 
factors, so that they cancel. Thus, we find that 



€SR> -{nH^)HAiNy 



(C7) 



is the consistency condition for linearized theory at second order. 
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APPENDIX D: SCALAR METRIC FLUCTUATIONS IN DE-SITTER 

In this section we comment that to second order in perturbation theory about pure de-Sitter (no-roll), 
the scalar sector is nontrivial. This can also be seen by e.g. examining the second order gauge fixing (45) 
in the main paper or the reduced equation of motion for Q, equation (48), which shows that the TT sector 
will mix with and source a nontrivial scalar sector at second order when esi? — > 0. It is perhaps worth 
emphasizing that the extra metric functions one can fix in flat, vacuum, spacetime and special spacetimes 
like de Sitter are generic, but obey field equations. Though they correspond to residual degrees of freedom 
specified only on an initial value surface, they are completely determined by equations of motion they also 
satisfy. Therefore any residual gauge fixing makes crucial use of the equations of motion to set additional 
terms to zero via a gauge transformation. 

To provide a formal alternate argument which however is far more compact and hopefully more convincing, 
consider the following. To linear order in perturbation theory about pure vacuum de-Sitter it is relatively 
straightforward that one can gauge away the scalar sector. Let us go through the proof and then broadly 
identify the impediment which would provent extending this procedure to second order. 

One may compactly prove this by writing by linearizing the linearized Einstein equations about 

ds^ = gabdx^dx'' = -dt'^ + cosh(i)2 [dx^ + sm'^{x)dn'^] , (Dl) 

which yields 

° Sgab - ^(jab ^ Sg- 26gab - gabSg = (D2) 
provided we choose the transverse Lorentz gauge, i.e. 



Sgab - ^Sg = 0. (D3) 



One then asks if one can impose tracelessness as well, and it turns out that making the additional choice 
Sg'i^i^ = 6gab + ^VaVh(5g is both consistent with the imposition of the Lorentz gauge (D3) and also satisfies 
the field equations (D2). However, the imposition of the Lorentz gauge along with tracelessness still leaves 
a residual gauge freedom in the form of gauge transformations which are harmonic functions of the operator 

□ +A, with which we can hope to eliminate 6gQ^. Not only does it turn that such transformations are 
possible, but the field equations expressed in this new gauge (combined with the maximal symmetry of de 
Sitter) also demand that (S^oo — 0- the complete, closed, slicing (Dl) the fact that the lapse fiuctuation 
SgQQ is zero corresponds to the fact that the eigenfunctions of the Laplacian on S'^ must be periodic (or, in 
other words, their eigenspectrum must be discrete). In other words, a positivity argument shows that the 
linear lapse fluctuation can also be gauged away. Thus tracelessness, transverseness, and the vanishing of 
the shift and lapse leaves two degrees of freedom, which are the two polarizations of the graviton. 

Therefore, at linear order in perturbation theory, the scalar sector is trivial (as is the vector-the proof is 
considerably simpler) modulo a linear gauge transformation and only the TT sector remains, which to this 
order is gauge invariant anyway. 

At second order, however, additional gravitational wave terms in TT-TT combinations (which transform 
as scalars) will act as sources for the scalar perturbations. These additional terms will enter as both negative 
definite and positive definite contributions because of the form of the perturbed Christoffel symbols (whose 
quadratic products and derivatives comprise the field equations), and this variable sign contribution is 
sufficient to in general get a nontrivial solution for the second order lapse. Therefore it is not possible in 
general to gauge it away at second order if one also insists on gauging away the second order shifts, so the 
second order scalar sector is not a gauge artifact. 
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